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The angular power spectrum of the cosmic microwave background (CMB) contains information 
on virtually all cosmological parameters of interest, including the geometry of the Universe (Q), 
the baryon density, the Hubble constant (h), the cosmological constant (A), the number of light 
neutrinos, the ionization history, and the amplitudes and spectral indices of the primordial scalar 
and tensor perturbation spectra. We review the imprint of each parameter on the CMB. Assuming 
only that the primordial perturbations were adiabatic, we use a covariance-matrix approach to 
estimate the precision with which these parameters can be determined by a CMB temperature map 
as a function of the fraction of sky mapped, the level of pixel noise, and the angular resolution. 
For example, with no prior information about any of the cosmological parameters, a full-sky CMB 
map with 0.5° angular resolution and a noise level of 15 per pixel can determine SI, h, and 
A with standard errors of ±0.1 or better, and provide determinations of other parameters which 
are inaccessible with traditional observations. Smaller beam sizes or prior information on some of 
the other parameters from other observations improves the sensitivity. The dependence on the the 
underlying cosmological model is discussed. 

98.70.V, 98.80.C 



I. INTRODUCTION 

One of the fundamental goals of observational cosmology today is measurement of the classical cosmological param- 
eters: the total density (or equivalently, the geometry) of the Universe, il; the cosmological constant A; the baryon 
density Qb] and the Hubble constant Hq. Accurate measurement of these quantities will test the cornerstones of 
the hot big-bang theory and will provide answers to some of the outstanding questions in cosmology. For example, 
determination of the geometry of the Universe will tell us the ultimate fate of the Universe and test the inflationary 
paradigm, while an independent check of fib can confirm the predictions of big-bang nucleosynthesis. 

In addition, parameters describing primordial perturbations are related to the origin of large-scale structure in 
the Universe and may shed light on a possible inflationary epoch. Perhaps the most important of these are the 
normalization Qs and spectral index ns of the primordial spectrum of scalar perturbations that gave rise to the 
observed structure. Inflation may produce a spectrum of gravity waves, quantified by an amplitude Qt and spectral 
index n-p. A neutrino species with a mass greater than an eV affects structure formation, so the number N^, of light 
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(meaning rriu ^ 1 eV) neutrinos is another cosmological parameter of importance. The ionization history of the 
Universe is also certainly related to the evolution of structure in the Universe. 

In this paper, we estimate how well cosmological parameters can be determined from a CMB temperature map. 
Since the initial detection of temperature anisotropics in the cosmic microwave background (CMB) by the COBE 
satellite 0, over a dozen other balloon-borne and ground-based experiments have announced anisotropy detections 
on smaller angular scales 1^. With the existence of anisotropics now firmly established, sights are shifting to accurate 
determination of the CMB power spectrum over a wide range of angular scales. Several technological advances, 
including improved amplifiers, interferometry, and long-duration balloon flights, hold great promise for high-precision 
measurements. Ultimately, a satellite with sub-degree angular resolution will provide a detailed map of the entire 
microwave sky in multiple frequency bands . 

A detailed map of the cosmic microwave background can potentially provide a wealth of information on the values 
of cosmological parameters. Roughly speaking, the amount of information in a map is proportional to the number 
of pixels on the sky, and this is inversely proportional to the square of the beam width. Thus, a map with a beam 
width of 0.5° will contain over 100 times as much information as COBE, which had a beam width of order 7°, and an 
0.1 "-resolution experiment would have, roughly speaking, lO'' times as much information! It should be no surprise, 
therefore, that a map with good angular resolution should be able to determine many more cosmological parameters 
than COBE, which really only constrains the normalization of the CMB power spectrum and the effective CMB 
spectral index at large angular scales. 

We consider an experiment which maps a given fraction of the sky with a given angular resolution and a given level 
of pixel noise. We use a covariance-matrix approach to evaluate the standard errors which would arise by fitting the 
power spectrum obtained in this experiment to all the unknown cosmological parameters. We display results for a 
range of realistic values for the fraction of sky covered, level of pixel noise, and angular resolution. Our results are 
quite promising: With minimal assumptions, realistic satellite experiments could potentially determine fi. A, and the 
inflationary observables to far greater precision than any traditional measurements. Furthermore, the information 
provided on other parameters will be competitive with (and with additional reasonable assumptions, superior to) 
current probes. Although we focus here only on models with primordial adiabatic perturbations, we are confident 
that if the perturbations turn out to be isocurvature, it will be evident in the temperature maps (and perhaps also 
in polarization maps, spectral distortions, and non-gaussian temperature distributions), and that similar results on 
parameter determination will apply. Indeed, recent calculations of the CMB power spectrum in defect models Q and 
in isocurvature models suggest that such models should be clearly distinguishable from the adiabatic case. Although 
we have satellite mapping experiments in mind, our results can also be applied to ground or balloon experiments, or 
to the combined results of several complementary measurements. 

An important issue facing any likelihood analysis is the choice of the space of models considered. Here we consider 
models with primordial adiabatic perturbations. Our space of models allows a cosmological constant, an open (or 
closed) Universe, tensor modes (with a free spectral index), variations in the baryon density and Hubble constant, 
tilted primordial spectra, and primordial spectra that deviate from pure power laws. We assume that the dark matter 
is cold; however, the CMB power spectrum is only slightly altered in mixed and hot dark-matter models and 
we allow the number of massless neutrinos to vary. Therefore, our conclusions on parameter determination will be 
virtually independent of the fraction of hot dark matter. 

In the following Section, we describe our calculation of the power spectrum. In Section III, we illustrate the effect 
of each cosmological parameter that we consider on the CMB spectrum. In Section IV, we discuss the covariance 
matrix. To illustrate, in Section V, we present results for the standard errors to the parameters that would be obtained 
assuming the true cosmological model is standard CDM. We also discuss how these results change if the underlying 
model differs from the canonical standard-CDM model. In Section VI, we discuss the validity of the covariance- 
matrix approach to the analysis. In Section VII, we make some concluding remarks and discuss some future areas of 
investigation. 

II. CALCULATION OF THE CMB SPECTRUM 

In many areas of astrophysics, it is difficult to make detailed quantitative predictions as properties of complex 
systems depend on non-linear physics of poorly measured and poorly understood phenomena. Fortunately, the early 
Universe was very simple and nearly uniform. The density fluctuations are all in the linear regime {5p/ p ^ 10^'*) 
and non-linear effects are unimportant. Different groups using different gauge choices and numerical algorithms make 
very similar predictions for CMB fluctuations for a given model. This simple linearity makes possible the detailed 
parameter determination that we describe in this paper. 

The CMB angular power spectrum C{0) is defined as 
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Ci9)^(^ — {rh) — {n)j, m.n = cos0, (1) 

where the angle brackets represent an ensemble average over all angles and observer positions. Here AT(n)/To is 
the fractional temperature fluctuation in the direction n, and the mean CMB temperature is Tq — 2.726 ± 0.010 K 
0. This power spectrum is conveniently expressed in terms of its multipolc moments Ci, defined by expanding the 
angular dependence in Legendre polynomials, Pe{x): 

C(0)=^-^QF,(cos0). (2) 

Given a model for structure formation, calculation of the multipole moments is straightforward and is accomplished by 
solution of the coupled system of Boltzmann equations for each particle species (i.e., photons, baryons, massless and 
possibly massive neutrinos, and cold dark matter) and Einstein equations for the evolution of the metric perturbations. 
The £ = 1 term is indistinguishable from the Doppler shift due to proper motion with respect to the microwave 
background rest frame and is conventionally ignored. For theories with gaussian initial perturbations, the set of C'i 
completely specifies the statistical properties of the theory. Since we can only observe from a single vantage point in the 
Universe, the observed multipole moments C°^^ will be distributed about the mean value Ci with a "cosmic variance" 

~ a/2/(2£ + 1)C^; no measurement can defeat this variance. Power-spectrum predictions and measurements are 
traditionally plotted as £{£ + l)Ce versus £. 

For the purposes of covariance-matrix evaluation, as well as for likelihood maximization |^,^ and Monte Carlo 
analysis, it is useful to have an algorithm for rapid evaluation of the CMB spectrum for a given set of cosmological 
parameters. We begin with a semi-analytic solution of the coupled Boltzmann, fluid, and Einstein equations developed 
by Hu and Sugiyama pO[ | for flat cold-dark-matter models, which we generalize to accommodate an open Universe, a 
cosmological constant, tensor modes, and reionization. The code is fast enough to enable likelihood analyses requiring 
tens of thousands of power-spectrum evaluations. We have checked that our semi- analytic calculation agrees with the 
results of a publicly available numerical code [ pl| for several parameters. Here we briefly describe our calculation. 

The multipole moments are expressed as 

Ci = Cf -I- Cf , (3) 

where is the contribution from scalar perturbations and Cj is the contribution from tensor modes. The scalar 
contribution is given by 

2 ^oo 

Cf = -/ dke\eiirjo,k)\^, (4) 
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where 770 is the conformal time today (the conformal time 77 = J dt/a with a the scale factor of the Universe normalized 
to unity at matter-radiation equality). The contribution of wavenumber k to the £th multipole moment is [|l0| 

rno 

6^(770, fc) - [6o + ^'](fc, J?* )ji(fc?7o-fc77*) + 0i(fc,??*)ji(fc??o-fcr7,)+ / di] ^ ^]ji{kriQ ~ ki]), (5) 

where Qq and Qi are the monopole and dipole perturbations of the photon distribution function, <I> and ^ are 
gravitational-potential perturbations in the Newtonian gauge, ji are spherical Bessel functions and j; their flrst 
derivatives, and a dot denotes derivative with respect to conformal time. Here 77* is the conformal time at decoupling. 
(See Ref. for more details.) The third term in this expression gives the integrated Sachs-Wolfe (ISW) effect: 
anisotropics are generated by time variations in the gravitational potentials along the line-of-sight path. Analytic fits 
to the gravitational potentials are given in Ref. jl^] , as are WKB solutions for the photon distributions in the tight- 
coupling regime, 0o and 0i. At decoupling, photon diffusion (Silk damping) damps photon perturbations on small 
angular scales [Q; the perturbations to the photon distribution functions are given by [9o-|-^'](77*) — [6o+4'](?7*)2?(fc), 
where the mean damping factor is 

V{k) = / fe-^^'^^^'^^^ dr]. (6) 
Jo 

Here f = XeneO'TCi/0.0 is the differential optical depth for Thomson scattering, Ue is the electron density, Xe is the 
ionization fraction, and ar is the Thomson cross section. The visibility function — the combination fe""^ — is the 
probability density that a photon last scattered at given conformal time, and is sharply peaked near the surface of 
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last scatter; semi-analytic fits are given in Ref. [T^]. As pointed out in Ref. [Q, photon polarization must be included 
to obtain the proper Silk-damping scale; the result is 

where 

is the scale factor normalized to 3/4 at baryon-radiation equality, with the fraction of critical density in baryons, 
Oq the fraction of critical density in non-relativistic matter (baryons and cold dark matter), and the fraction of the 
total radiation density contributed by massless neutrinos. Our numerical evaluation of these expressions reproduces 
the power spectrum obtained from Boltzmann codes to an accuracy of a few percent for standard CDM. 

Analytic approximations to the CMB anisotropy due to tensor modes (gravity waves) are given in Refs. 
The contribution to each multipole moment of the CMB power spectrum is 

Cl = S6n^T7^ r dkPTik)\Feik)\\ (9) 
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where Pt oc fc"^^^ is the initial power spectrum of tensor perturbations and Fg is given by 

with fccq defined as the wavenumber of the mode which enters the horizon at matter-radiation equality. The fitting 
function w{rj) describes the evolution of the gravity- wave mode function through the transition between the radiation- 
dominated and matter-dominated epochs, and T{k, rj) is a transfer function describing the evolution of the tensor-mode 
amplitude. Good analytic fits to these two functions are given by 



w(r/) = exp (-0.27]°-5'"^) (11) 

2 r 

T(y,77) = ^ e-4^'(l + 1.342/ + 2.5y2)i/2 + i_g-V . (12) 
a L J 

These approximations match numerical results to one percent well past £ — 100, where the tensor contribution to the 
multipoles drops to a small fraction of the scalar contribution. 

Eqs. (^ and (H) are difficult to evaluate numerically because of the oscillatory spherical Bessel functions in the 
integrand. Asymptotic expansions, a Bessel-function cache, and various interpolation techniques further speed eval- 
uation of the integrals. We calculate every 40th multipole (more for £ < 100) and perform a cubic spline to recover 
the entire spectrum. 

We consider models which are well-described by a power-law spectrum of metric perturbations over the range of 
scales affecting CMB anisotropics. This class includes all inflation models. For the scalar perturbations, we also allow 
a deviation from power-law behavior and parameterize the power spectrum as p6| 

\ ns+a \n(k/ks) 



P{k) « j , (13) 

where ks is the normalization scale at which the power law index ns is defined. The parameter a quantifies the 
deviation from the power law, or the "running" of the spectral index. Realistic inflation models can produce values 
of a large enough to change the multipole moments by as much as 5%. For the tensor spectrum, we assume a 
pure power-law spectrum with spectral index nx. In principle, ut can run with scale as well, but because of the 
comparatively small amount of information contained in the tensor multipole moments, the CMB constraint on the 
index riT is weak, and the running-index effect for the tensor perturbations is negligible. 

Extensions of this basic cosmological model are incorporated through various fitting formulas. In a cosmological- 
constant (A) Universe, the gravitational potential $ begins to vary at low redshift when the Universe becomes 
cosmological-constant dominated, and this leads to a contribution to the anisotropy at large angles from the integrated 
Sachs- Wolfe (ISW) effect. In a flat Universe (that is, fig + A = 1, where A is the cosmological constant in units of 
critical density), this is approximated by multiplying the multipole moments by a factor [1 + g{h)/£] ||l^,f8|, where 
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is the time dependence of the potential, and H = a/a \s the Hubble parameter. This approximation slightly over- 
estimates the lowest few multipole moments, but this large-angle ISW effect is generally not a large fraction of the 
total anisotropy, and the lowest multipole moments have a limited statistical significance. For A < 0.7, (/(A) can be 
approximated by 



A 



5(A) 0.637 (^^-^ j . (16) 

An additional effect of a cosmological constant is a shift in the conformal distance to the surface of last scatter (even 
with the mass density Vl^h'^ held fixed), which we account for by multiplying the current conformal time 770 by the 
correction factor 1 + 0.085 ln(l - A) [|l9|. 

Generalization to an open Universe is somewhat more complicated because several different effects contribute to 
the anisotropy [Q. The angular scale subtended by the horizon at the surface of last scatter scales as VL^/^ where 

= ilo + A is the total density (in units of critical density) of the Universe [|2^ . Therefore, the multipole moments in an 
open Universe are related to those in a flat Universe approximately by Di{^) ~ £'^qi/2(S1 — 1) with = £{£+ l)Ce. 
In other words, the CMB spectrum in an open Universe resembles that in a flat Universe with the same matter 
density, but shifted to smaller angular scales. A large-angle ISW effect arises from the evolution of the gravitational 
potentials, although the function g{^l) differs from that in a cosmological-constant Universe [|l8|. In addition, the 
lowest multipole moments probe scales comparable to or larger than the curvature scale, so these moments are 
suppressed, due heuristically to the exponential growth of volume in an open Universe at large distances. Finally, 
some ambiguity exists as to the correct generalization of a power-law spectrum to an open Universe. Naive power laws 
of volume, wavenumber, or eigenvalue of the Laplace operator differ in an open Universe |p^ , as do spectra predicted 
by various open- Universe inflationary scenarios [ pl| . However, these power laws differ only in their predictions for 
the lowest multipole moments, which have little statistical weight; for definiteness, we use the predictions of a specific 
inflationary scenario j22|. A good fit to these effects (for Q > 0.1) is provided by multiplying the multipole moments 
by 



-0.3 f/^curv g(^) 



£+1/2' 



where £curv = tt-^/ (1 — n)/il is the multipole corresponding to the curvature scale of the Universe, and 



0.817 



(17) 



g(r!) ^ 4.5 (^-^ j , (18) 
for n > 0.1. 

If the Universe has experienced significant reionization between recombination and today, then a fraction 1 — e^'^rcio„ 
of the CMB photons have scattered since recombination, where Treion is the optical depth to the epoch of recombination. 
If the Universe becomes reionized at a redshift Zreion with a constant ionization fraction x^, then the optical depth is 
Treion — 0.04:Qbh^~^^^Xe[{l + ^roion)^^^ " 1], where h IS the Hubble parameter in units of 100 km sec~^ Mpc~^. The 
precise effects of reionization depend on the baryon density, Hubble parameter, and the ionization history. However, 
as illustrated in Ref. (see Fig. 3 therein) , the effects of reionization are fairly accurately quantified solely in terms 
of Ticion- Compton scattering is an isotropizing process, so the multipole moments on angular scales smaller than 
those subtended by the horizon at the epoch of reionization are suppressed by a factor e~'^'^"''°'^ , while those on larger 
angular scales are unaffected. We interpolate between the asymptotic effects of reionization on small and large angular 
scales by multiplying the multipole moments by 



exp 



reion 



1 + {(riicion/voV 



(19) 



where rj^eion is the conformal time at reionization. In addition, reionization also induces a broad Doppler peak 
centered near £ ~ ryo /?/roion | |2^ , |2^ , but this secondary peak is shallow and we do not include it in the power-spectrum 
calculation. 
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Between the surface of last scatter and the present, several other physical processes, besides reionization, produce 
new CMB fluctuations and smear out primordial fluctuations |^5| : gravitational lensing lowers the amplitude of the 
spectral peaks and fills in the valleys in the spectrum |2^; the non-linear growth of structure produces additional 
small-scale fluctuations j27j; the scattering of photons off of hot gas in clusters and superclusters produces both 
thermal and non-thermal cosmic microwave fluctuations ; and second-order effects in a reionized Universe also 

produce additional small-scale fluctuations [^,0 . These non-linear effects are relatively small and typically produce 
only ©(/iK) changes in the microwave multipolcs. However, they are systematic. If they are not included in an 
analysis of a full-sky CMB map, they will lead to systematic errors in parameter estimation. We do not include these 
effects in our sensitivity analyses as they are unlikely to alter the size and shape of the error ellipsoid. It will be 
important to include these effects in any analysis of a future all-sky CMB map. 

III. COSMOLOGICAL PARAMETERS AND THE CMB SPECTRUM 

The suite of cosmological models that we consider all make broadly similar predictions for the CMB spectrum: 
the fluctuations on large angular scales are nearly scale-invariant and are primarily due to large-scale variations in 
the gravitational potential at the surface of last scatter, while on small scales the fluctuations are primarily due to 
variations in the velocity and density of the baryon-photon fluid at the surface of last scatter. The details of the 
spectrum, however, depend sensitively on properties of the Universe: its geometry, its size, the baryon density, the 
matter density, and the shape of the primordial fluctuation spectrum. In this Section, we discuss each parameter that 
we consider and illustrate its most salient effect on the CMB spectrum. Fig. |^ illustrates the following discussion. 

The first Doppler peak occurs at the angular scale subtended by the sonic horizon at the surface of last scatter. 
Since the photon energy density exceeds the baryon energy density at that epoch, the sound speed of the Universe is 
close to c/\/3, so that the sonic horizon corresponds to a nearly fixed physical scale. The angular scale subtended by 
this fixed physical scale will depend on the geometry of the Universe. In an open Universe, the angular scale subtended 
by an object of fixed diameter at fixed large redshift scales as O. On the other hand, the causal horizon at last scatter 
is actually times as large in an open Universe as it is in a flat Universe. Thus, to a first approximation, the 

flat-Universe CMB spectrum is stretched by a factor fi^/^ to smaller angular scales in an open Universe. 

Increasing the baryon density, flbh^, reduces the pressure at the surface of last scatter and therefore increases the 
anisotropy at the surface of last scatter. This reduction in pressure also lowers the sound speed of the baryon-photon 
fluid, which alters the location and spacing of the Doppler peaks. Increasing the matter density, fio/i^, shifts matter- 
radiation equality to a higher redshift. This reduces the early-ISW contribution to the spectrum and lowers and 
narrows the flrst Doppler peak. If we knew that A = 0, then the combination of these three effects (pressure, sound 
speed, and redshift of matter-radiation equality) would be sufflcient to enable a determination of f2oj ^6, and h from 
the CMB spectrum. 

The cosmological constant introduces a near degeneracy in parameter determination. Bond ct al. plj stressed that 
the CMB spectrum changed little if A was varied while Uoh^ and fl^h^ were held fixed in a fiat Universe. Changing 
A, however, does alter the size of the Universe. The conformal distance from the present back to the surface of 
last scatter is smaller in a A-dominated flat Universe than in a matter-dominated flat Universe. Thus, increasing A 
shifts the Doppler peak to larger angular scales, the opposite effect of lower ^Iq. This effect, along with the late-time 
ISW effect induced by A, breaks the degeneracy and enables an independent determination of all of the cosmological 
parameters directly from an all-sky high-resolution CMB map. 

The value of iV^, the effective number of noninteracting relativistic degrees of freedom (in standard CDM, this is 
equal to three for the three light-neutrino species), also shifts the epoch of matter-radiation equality and thus the 
height of the first Doppler peak as discussed above. In addition, if N^, is changed, the value of the anisotropic stress 
at early times — before the Universe is fully matter dominated — is altered, and this has a slight effect on the ISW 
contribution to the rise of the first Doppler peak. 

The tensor-mode contribution to the multipole moments simply adds in quadrature with the scalar-mode contri- 
bution since there is no phase correlation between them. The amplitude of the tensor modes is parameterized by 
r — Q^/Qg, the ratio of the squares of the tensor and scalar contributions to the quadrupole moment. [] The index 
riT is defined so that the tensor-mode spectrum is roughly fiat at large angular scales for ht = 0; it falls steeply near 
the rise to the first Doppler peak. Thus, tensor modes may contribute to the anisotropy at large scales, but they will 



^Note that this definition differs from that in Ref. 
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have little or no effect on the structure of the Doppler peaks. Increasing the tensor spectral index, ut, increases the 
contribution at small angular scales relative to those at larger angles. 

The overall normalization, Q, raises or lowers the spectrum uniformly. The effect of the scalar spectral index is 
similarly simple: if ns is increased there is more power on small scales and vice versa. The effects of a are obvious 
from Eq. (13). Finally, the effects of reionization have been discussed in the previous Section. 



IV. ERROR ESTIMATES 



We consider an experiment which maps a fraction /sky of the sky with a gaussian beam with full width at half 
maximum 6'fwhm and a pixel noise CTpix = s/V^pix, where s is the detector sensitivity and tpix is the time spent 
observing each 0fwhm x ^fwhm pixel. We adopt the inverse weight per solid angle, = (crpix6'fwhm/'7o)^, as a measure 
of noise that is pixel-size independent . Current state-of-the-art detectors achieve sensitivities of s = 200 /iK v^sec, 
corresponding to an inverse weight ofw^^ ~ 2 x 10~^^ for a one- year experiment. Realistically, however, foregrounds 
and other systematic effects may increase the effective noise level; conservatively, will likely fall in the range 
(0.9 - 4) X 10"". Treating the pixel noise as gaussian and ignoring any correlations between pixels, estimates of Cg 
can be approximated as normal distributions with a standard error 



2 


1/2 


[(2£+l)/skyJ 





(20) 



where ai, = 7.42 x 10 (^'fwhm/l°)-0 Note that Eq. ^ applies only if the entire sky has been mapped and then 
a fraction 1 — /sky has been subtracted. On the other hand, if only a fraction /sky of the sky is mapped, then the 
integration time per pixel increases by a factor of f^y, and w^^ should be replaced by iw^^/sky p^ . 

In Fig. ^, we show simulated data that might be obtained with a CMB mapping experiment, given an underlying 
cosmological model of "standard CDM" (see the following section). The "Cosmic Variance" panel illustrates the 
multipole moments that would be measured by an ideal experiment (i.e., perfect angular resolution and no pixel 
noise); the scatter is due only to cosmic variance. The top-right and bottom-left panels show multipole moments that 
might be measured by full-sky mapping experiments with a realistic level of pixel noise and angular resolutions of 
0.1° and 0.3°, respectively. The cosmic variance slightly increases the errors at lower i, while the finite beam width 



/0.3°) in the lower-left plot. The lower-right panel shows the 



is evident in the increased noise at (£/700) > (0fwhm 
moments from the lower-left panel after the total signal is smoothed with a gaussian window of width ^/20. This 
illustrates that although the individual moments may be quite noisy, an experiment with a beam width of 0.3° can 
still use the information in the location and shape of the third peak in parameter estimation. An experiment with this 
size beam can extract useful information out to f ^ 900, although it can not accurately measure the individual values 
of these high £ multipoles. The smoothing here is used for display and is not the optimal approach for parameter 
estimation. 

We now wish to determine the precision with which a given CMB temperature map will be able to determine the 
various cosmological parameters. The answer to this question will depend not only on the experimental arrangement, 
but also on the correct underlying cosmological parameters which we seek to determine. For any given set of cos- 
mological parameters, s = {fl,ni,h^,h,A,ns,r,nT,cx,TreiomQ,N^}, the multipole moments, Ci{s), can be calculated 
as described above. Suppose that the true parameters which describe the Universe are Sq. If the probability for 
observing each multipole moment, C^^^, is nearly a gaussian centered at Ci with standard error ag, and 0fwhm <C 1 so 
that the largest multipole moments sampled are £ 3> L then the probability distribution for observing a CMB power 
spectrum which is best fit by the parameters s is |3|,M,^ 



P(s) cx exp 



1 



(s - So) • [a] ■ (s - So) 



where the curvature matrix [a] is given approximately by 

1 rac,(so)ac,(so) 



E 



ds-i 



(21) 



(22) 



We thank a referee for pointing out an error in this equation in an earlier draft. 
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As discussed in Ref. |g2[, the covariance matrix [C] = [a] ^ gives an estimate of the standard errors that would 
be obtained from a maximum-Ukehhood fit to data: the standard error in measuring the parameter Si (obtained 

1 /2 

by integrating over all the other parameters) is approximately CJ . Prior information about the values of some of 
the parameters — from other observations or by assumption — is easily included. In the simplest case, if some of the 
parameters are known, then the covariance matrix for the others is determined by inverting the submatrix of the 
undetermined parameters. For example, if all parameters are fixed except for si, the standard error in Si is simply 

-1/2 

Previous authors have investigated the sensitivity of a given experimental configuration to some small subset of 
the parameters we investigate here. For example, Knox investigated the sensitivity of mapping experiments to the 
inflationary parameters, ng, tit, and r, but assumed all other parameters (including fli, and h) were known 
Similarly, Hinshaw, Bennett, and Kogut investigated the sensitivity to fib assuming all other parameters were fixed 
p7| . These were Monte Carlo studies which mapped the peak of the likelihood function. Another technique is to 
repeatedly simulate an experimental measurement of a given underlying theory, maximize the likelihood in each case 
and see how well the underlying parameters are reproduced j^. Such calculations require numerous evaluations of 
the CMB spectrum, so the results have been limited to a small range of experimental configurations. If any of these 
analyses are limited to a small subset of cosmological parameters, they do not investigate the possible correlation 
with other undetermined parameters and will therefore overestimate the capability of the experiment to measure the 
parameters under consideration. 

The covariance-matrix approach has the advantage that numerous experimental configurations and correlations 
between all the unknown cosmological parameters can be investigated with minimal computational effort. For example, 
if there are A'^ undetermined parameters, then we need only A^ + 1 evaluations of the Cg's to calculate the partial 
derivatives in Eq. (21). Once these are evaluated, the curvature matrix for any combination of and ^fwhm 
for /sky = 1 can be obtained trivially. The results are generalized to /sky < 1 by multiplying the results for the 
curvature matrix by f~^y [c.f., Eqs. (^0|) and (p2[)]. Furthermore, the covariance matrix includes all correlations 
between parameters. Therefore, our results reproduce and generalize those in Refs. [^,^,^, and we comment on this 
further below. 



V. COVARIANCE-MATRIX RESULTS 



As discussed above, the sensitivity of a CMB map to cosmological parameters will depend not only on the experi- 
ment, but also on the underlying parameters themselves. For illustration, we show results for a range of experimental 
parameters under the assumption that the underlying cosmological parameters take on the "standard-CDM" values. 
So = {1, 0.01, 0.5, 0, 1, 0, 0, 0, 0, QcoBE, 3}, where Qcobe = 20 is the COBE normalization Q. (It assumes a 
Harrison-Zeldovich primordial spectrum, no tensor modes, no cosmological constant, a flat Universe, and the central 
big-bang nucleosynthesis value for the baryon-to-photon ratio.) After presenting results for this assumed cosmological 
model, we briefly discuss how the results will be altered for different cosmological models. 

With the eleven undetermined cosmological parameters we survey here — some of which are better determined by 
experiment than others — there are an endless number of combinations that could conceivably be investigated. Instead 
of running through all possible permutations, we present results for the standard errors that can be obtained with two 
extreme sets of assumptions. First, we consider the case where none of the parameters are known. Then we consider 
the results under the most optimistic assumption that all of the other parameters, except the normalization (which 
will never be determined more accurately by any other observations), are fixed. Realistically, prior information on 
some of the parameters will be available, so the standard errors will fall between these two extremes. 

Figs. ^ and |^ show the standard errors for various parameters that can be obtained with a full-sky mapping 
experiment as a function of the beam width 0fwhm for noise levels — 2 x 10^^'°^, 9 x 10~^^, and 4 x 10" (from 
lower to upper curves). The underlying model is "standard CDM." The solid curves are the sensitivities attainable 
with no prior assumptions about the values of any of the other cosmological parameters. The dotted curves are the 
sensitivities that would be attainable assuming that all other cosmological parameters, except the normalization (Q), 
were fixed. The analogous results for a mapping experiment which covers only a fraction /sky of the sky can be 

obtained by scahng by /sky^^ [c-f-, Eq. (f^)]. 



A. The Total Density and Cosmological Constant 

The results for O were discussed in Ref. [Q. From the O panel in Fig. ^, it should be clear that a CMB mapping 
experiment with sub-degree resolution could potentially determine fl to better than 10% with minimal assumptions. 
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and perhaps better than 1% with prior information on other cosmological parameters. This would be far more precise 
than any conventional measurement of fl. Furthermore, unlike mass inventories which measure only the matter density 
fio, this measurement includes the contribution to the density from a cosmological constant (i.e., vacuum energy) and 
therefore directly probes the geometry of the Universe. This determination follows from the angular location of the 
first Dopplcr peak. Therefore, our results show that if the Doppler peak is found to be at I' ~ 200, it will suggest a 
value of 17 = 1 to within a few percent of unity. This result will be independent of the values of other cosmological 
parameters and will therefore be the most precise test for the flatness of the Universe and thus a direct test of the 
inflationary hypothesis. Numerical calculations suggest that the effect of geometry on the CMB spectrum may be 
slightly more dramatic than indicated by our semi-analytic algorithm.]^ If so, our final results on the sensitivity to 
are a conservative estimate. 

The sensitivity to A is similar. Currently, the strongest bounds to the cosmological constant come from gravitational- 



lensing statistics |39 which only constrain A to be less than 0.5. Measurement of the deceleration parameter, go = 
r2o/2 — A, could provide some information on A, but the measurements are tricky, and the result will depend on the 
matter density. On the other hand, a CMB mapping experiment should provide a measurement of Lambda to better 
than ±0.1, which will easily distinguish between a A-dominated Universe and either an open or flat matter-dominated 
Universe. 

B. The Baryon Density and Hubble Parameter 

The current range for the baryon-to-photon ratio allowed by big-bang nucleosynthesis (BBN) is 0.0075 ^ ilt,h^ < 
0.024 1^^. This gives Qb ^ 0.1 for the range of acceptable values of h, which implies that if = 1, as suggested by 
inflationary theory (or even if 0.3 as suggested by cluster dynamics), then the bulk of the mass in the Universe 
must be nonbaryonic. On the other hand, x-ray-cluster measurements might be suggesting that the observed baryon 
density is too high to be consistent with BBN H] ; this becomes especially intriguing given the recent measurement of 
a large primordial deuterium abundance in quasar absorption spectra ||42f| . The range in the BBN prediction can be 
traced primarily to uncertainties in the primordial elemental abundances. There is, of course, also some question as 
to whether the x-ray-cluster measurements actually probe the universal baryon density. Clearly, it would be desirable 
to have an independent measurement of flhh^. The flhh^ panel in Fig. |^ shows that the CMB should provide such 
complementary information. The implications of CMB maps for the baryon density depend quite sensitively on the 
experiment. As long as 6'fwhm ^ 0.5, the CMB should (with minimal assumptions) at least be able to rule out a 
baryon-dominated Universe {fib ^ 0.3) and therefore confirm the predictions of BBN. With angular resolutions that 
approach 0.1° (which might be achievable, for example, with a ground-based interferometry map p3[ to complement 
a satellite map), a CMB map would provide limits to the baryon-to-photon ratio that were competitive with BBN. 
Furthermore, if other parameters can be fixed, the CMB might be able to restrict ri^/i^ to a small fraction of the 
range currently allowed by BBN. 

Current state-of-the-art measurements of the Hubble parameter approach precisions of roughly 10%, and due 
to systematic uncertainties in the distance ladder, it is unlikely that any determinations in the foreseeable future 
will be able to improve upon this result. The panel for h in Fig. || shows that, even with minimal assumptions, a 
mapping experiment with angular resolution better than 0.5° will provide a competitive measurement; with additional 
assumptions, a much more precise determination is possible. It should also be noted that the CMB provides a 
measurement of the Hubble parameter which is entirely independent of the distance ladder or any cosmological 
distance determination. 

As a technical aside, we mention that in calculating the curvature matrix, Eq. (p2[), we choose Qoh^ as an inde- 
pendent parameter instead of h, and then transform the curvature matrix back to the displayed parameters. The 
reason for this choice is that the power spectrum depends on h only indirectly through the quantities Oq/i^ and flbh^, 
and the linear approximation to the change in the spectrum in Eq. ( p2[ ) is more accurate for the parameter floh^. 
This parameter choice also explicitly accounts for the approximate degeneracy between models with the same value 
of Qoh^ but differing A [|l[. 



^We thank a referee for pointing this out. 
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C. Reionization 



As discussed above, the effects of reionization can be quantified, to a first approximation, by Treion, the optical depth 
to the surface of last scatter, and there are several arguments which suggest Ti-oion ^ 1 First of all, significant 

reionization would lead to anisotropies on arcminute scales due to the Vishniac effect |2^, or to spectral (Compton-y) 
distortions of the CMB |Q . Order-of-magnitude estimates for the values of Treion expected in adiabatic models based 
on Press-Schechter estimates of the fraction of mass in collapsed objects as a function of redshift suggest that Treion 
is probably less than unity [ pOp^ ]. Moreover, the numerous detections of anisotropy at the degree scale Q also show 
an absence of excessive reionization. Assuming complete reionization at a redshift Zioion, the optical depth with our 
standard-CDM values is Treion — 0.001 ^^fe^n' '^roion ^ 1 corresponds to Zrcion ^ 100. 

The Treion panel of Fig. ^ illustrates that, with minimal assumptions, any map with sub-degree angular resolution 
will probe the ionization history (i.e., Zroion ^ 1000), and maps with resolutions better than a half degree can restrict 
the optical depth to 0.5 or less. While different ionization histories with the same total optical depth can give different 
power spectra, as long as the reionization is not too severe, simple damping of the primary anisotropies is always 
the dominant effect. The lower curves, assuming other parameters are fixed except for Q, are flat because at scales 
smaller than 2°, the effects of Treion are precisely degenerate with a shift in Q. The lower curves nearly coincide for 
the different noise levels because all of the leverage in distinguishing Treion comes at low £ where the degeneracy with 
Q is broken, and at these scales the cosmic variance dominates the measurement errors. 

Although temperature maps alone may not provide a stringent probe of the ionization history, polarization maps 
may provide additional constraints The polarization produced at recombination is generally small, but that 
produced during reionization can be much larger. Heuristically, the temperature anisotropy which is damped by 
reionization goes into polarization. Therefore, it is likely that polarization maps will be able to better constrain Treion 
when used in conjunction with temperature maps. 



D. Neutrinos 



We have also investigated the sensitivity of CMB anisotropies to N^, the effective number of neutrino degrees of 
freedom at decoupling. The number of non-interacting relativistic species affects the CMB spectrum by changing the 
time of matter-radiation equality, although this cannot be distinguished from the same effect due to changes in h, 
rio, and A. However, neutrinos (and other non- interacting degrees of freedom which are relativistic at decoupling) 
free stream and therefore have a unique effect on the growth of potential perturbations. This will be reflected in the 
detailed shape of the CMB spectrum, especially from the contribution of the early-time ISW effect. In standard CDM, 
there are the three light-neutrino species. However, some particle-physics models predict the existence of additional 
very light particles which would exist in abundance in the Universe. Furthermore, if one of the light neutrinos has a 
mass greater than an eV, as suggested by mixed dark-matter models p6| and possibly by the Los Alamos experiment 
, then it would be nonrelativistic at decoupling so the effective number of neutrinos measured by the CMB would 
be iVy < 3.^ These limits would be similar to limits on the number of relativistic noninteracting species from BBN. 
However, at the time of BBN, any particle with mass less than an MeV would be relativistic, whereas at decoupling, 
only those with masses less than an eV would be relativistic, so the quantities probed by BBN and by the CMB are 
somewhat different. 

The panel for in Fig. ^ shows the sensitivity of CMB anisotropies to variations in the effective number of non- 
interacting nonrelativistic species at decoupling. When one takes into account systematic uncertainties in primordial 
elemental abundances, BBN constrains the effective number of relativistic (i.e., less than a few MeV) neutrino species 
to be less than 3.9 ]40[ |. Fig. 1 illustrates that any mapping experiment with angular resolution better than 0.5° 
should provide complementary information; if other parameters can be determined or constrained, then the CMB has 
the potential to provide a much more precise probe of the number of light neutrinos at the decoupling epoch. 



*In such a case, the massive neutrino would have additional effects on the CMB Q|. Although we have not included these 
effects, our analysis still probes variations in A'^^, and our results are suggestive of the sensitivity of CMB anisotropies to a 
massive neutrino. 
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E. Inflationary Observables 



We have also studied the precision with which the inflationary observables, ns, nx, and r, can be probed. Inflation 
predicts relations between the scalar spectral index ns, the tensor spectral index ut, and the ratio r [Q. Therefore, 
precise measurement of these parameters provides a test of inflationary cosmology and perhaps probes the inflaton 
potential jES] . 

Knox [03| performed a Monte Carlo calculation to address the question of how accurately CMB anisotropics can 
measure the inflationary observables assuming all other cosmological parameters were known. Here, we generalize 
the results to a broader range of pixel noises and beam widths and take into account the uncertainties in all other 
cosmological parameters through the covariance matrix. 

In Fig. ^, we show the standard errors on the inflationary observables that could be obtained with mapping 
experiments with various levels of pixel noise and beam widths. The parameters of the underlying model used here 
are the same "standard-CDM" parameters used in Figs. || and U, except here we set r = Q^/Q^ = 0.28, ns = 0.94, 
and nx = —0.04. We do so for two reasons: First, the tensor spectral index is unconstrained without a tensor 
contribution; second, these parameters will facilitate comparison with the results of Ref. |33|. The solid curves are 
the standard errors that would be obtained with no assumptions about the values of these or any other of the 
cosmological parameters. The dotted curves are the standard errors that would be attainable by fitting to only these 
four inflationary observables and assuming all other cosmological parameters are known. (Note that this differs from 
the dotted curves in Figs. I and |.) 

The dotted curves in Fig. |^ with a beam width of 0.33° are in good numerical agreement with the results of 
Ref. [Q. This verifies that the covariance-matrix and Monte Carlo calculations agree. Next, note that unless the other 
cosmological parameters can be determined (or are fixed by assumption), the results of Ref. |]33|] for the sensitivities 
of CMB anisotropy maps to the inflationary observables are very optimistic. In particular, temperature maps will be 
unable to provide any useful constraint to r and n^ (and it will be impossible to reconstruct the inflaton potential) 
unless the other parameters can be measured independently. However, if the classical cosmological parameters can be 
determined by other means (or fixed by assumption), the dotted curves in Fig. ^ show that fairly precise information 
about the inflaton potential will be attainable. CMB polarization maps may provide another avenue towards improved 
determination of the inflationary observables ||,|5^ . 

The flatness of the dotted curves for r and nx in Fig. || is due to the fact that the contribution of the tensor modes 
to the CMB anisotropy drops rapidly on angular scales smaller than roughly a degree. The solid curves decrease with 
^fwhm because the other cosmological parameters (e.g., ftbh^ and h) become determined with much greater precision 
as the angular resolution is improved. 

Of course, the precision with which the normalization of the perturbation spectrum can be measured with CMB 
anisotropies (even current COBE measurements) is — and will continue to be — unrivaled by traditional cosmological 
observations. Galactic surveys probe only the distribution of visible mass, and the distribution of dark matter could be 
signiflcantly different (this is the notion of biasing). The dotted figures in the panel for Q in Fig. |^ are the sensitivities 
that would be obtained assuming all other parameters were known. This standard error would be slightly larger if 
there were no tensor modes included, because as r is increased (with the overall normalization Q held fixed), the 
scalar contribution is decreased. Therefore, tensor modes decrease the anisotropy on smaller angular scales, the signal 
to noise is smaller, and the sensitivity to Q (and other parameters) is slightly decreased. The effect of variations in 
other underlying-model parameters on our results is discussed further below. 



F. What If The Underlying Model Is Diff'erent? 

Now we consider what might be expected if the underlying theory differs from that assumed here. Generally, the 
parameter determination will be less precise in models in which there is less cosmological anisotropy, as refiected in 
Eq. (§). 

What happens if the normalization differs from the central COBE value we have adopted here? The normalization 
raises or lowers all the multipole moments; therefore, from Eq. (pO|), the effect of replacing Q with Q' is equivalent 
to replacing w"^ with w~'^{Q/Q'). In Figs. ^, ^, and the solid curves, which are spread over values of that 
differ by more than an order of magnitude, are all relatively close. On the other hand, the uncertainty in the COBE 
normalization is 0(10%). Therefore, our results are insensitive to the uncertainty in the normalization of the power 
spectrum. 

If there is a significant contribution to the COBE-scale anisotropy from tensor modes, then the normalization of the 
scalar power spectrum is lower, the Doppler peaks will be lower, and parameter determinations that depend on the 
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Doppler-peak structure will be diluted accordingly. On the other hand, the tensor spectral index, which is important 
for testing inflationary models, will be better determined. 

Similarly, reionization damps structure on Doppler-peak angular scales, so if there is a significant amount of reion- 
ization, then much of the information in the CMB will be obscured. On the other hand, there are several indications 
summarized above that damping due to reionization is not dramatic. In Ref. [ ^2| , we displayed (in Fig. 2 therein) 
results for the standard error in Q for a model with Tioion = 0.5. As expected, the standard error is larger (but by no 
more than a factor of two) than in a model with no reionization. 

If A is non-zero, h is small, or ili,h^ is large, then the signal-to-noise should increase and there will be more 
information in the CMB. If the scalar spectral index is ns > 1, then the Doppler peaks will be higher, but in the 
more likely case (that predicted by inflation) , ns will be slightly smaller than unity. This would slightly decrease the 
errors. 

If is less than unity, then the Doppler peak (and all the information encoded therein) is shifted to smaller angular 
scales. Thus one might expect parameter determinations to become less precise if < 1. By explicit numerical 
calculation, we find that for = 0.5 (with all other parameters given by the "standard-CDM" values), our estimates 
for the standard error for SI is at the same level as the estimate we obtained for = 1. Therefore, even if is as small 
as 0.3, our basic conclusions that Q can be determined to ±0.1 with minimal assumptions are valid. The sensitivity 
to some of the other parameters, in particular h, is (not surprisingly) significantly degraded in an open Universe. 



VI. GAUSSIAN APPROXIMATION TO THE LIKELIHOOD FUNCTION 



It is an implicit assumption of the covariance-matrix analysis that the likelihood function has an approximate 
gaussian form within a sufficiently large neighborhood of the maximum-likelihood point. Eq. (^ ) only approximates 
the likelihood function in a sufficiently small neighborhood around the maximum. The detailed functional form of the 
likelihood function is given in Ref. [ p3[ . If the likelihood function fails to be sufficiently gaussian near the maximum, 
then the covariance-matrix method is not guaranteed to produce an estimate for the standard error, and a more 
involved (Monte Carlo) analysis would be essential. Therefore, in the following we indicate the applicability of the 
gaussian assumption for the likelihood. 

First, we note that our parameterization has the property that the individual parameters are approximately inde- 
pendent. This is suggested by direct examination of the eigenvectors of the covariance matrix. This is also supported 
by preliminary Monte Carlo results. Therefore it is simplest to examine the behavior of the likelihood as a function of 
individual parameters in order to determine if a parabolic approximation to ln(£) (the log-likelihood) is admissible. 

In Fig. ^ we display the dependence of ln(£) on several parameters of interest. In this example, we used w^^ = 
9 X 10~^^ and 6'fwhm = 0.25. As is clear from this Figure, the functional forms are well fit by parabolic approximations, 
within regions of size ~ Str around the maximum point. This is sufficient to apply the covariance-matrix analysis to 
the determination of the standard errors, and our analysis above is justified. 

Although we have not done an exhaustive survey of the likelihood contours in the eleven-dimensional parameter 
space, Fig. |^ also suggests that there are no local maxima anywhere near the true maximum. Therefore, fitting 
routines will probably not be troubled by local maxima. This also suggests, then, that there will be no degeneracy 
between various cosmological models with a CMB map (in contrast to the conclusions of Ref. |31|), unless the models 
are dramatically different. In this event (which we consider unlikely), one would then be forced to choose between two 
quite different models, and it is probable that additional data would determine which of the two models is correct. 



VII. CONCLUSIONS 



We have used a covariance-matrix approach to estimate the precision with which eleven cosmological parameters 
of interest could be determined with a CMB temperature map. We used realistic estimates for the pixel noise and 
angular resolution and quantified the dependence on the assumptions about various cosmological parameters that 
would go into the analysis. The most interesting result is for Q: With only the minimal assumption of primordial 
adiabatic perturbations, proposed CMB satellite experiments |^ could potentially measure Q to 0(5%). With prior 
information on the values of other cosmological parameters possibly attainable in the forthcoming years, f2 might be 
determined to better than 1%. This would provide an entirely new and independent determination of fl and would 
be far more accurate than the values given by any traditional cosmological observations. Furthermore, typical mass 
inventories yield only the matter density. Therefore, they tell us nothing about the geometry of the Universe if the 
cosmological constant is nonzero. A generic prediction of infiation is a flat Universe; therefore, locating the Doppler 
peak will provide a crucial test of the inflationary hypothesis. 
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CMB temperature maps will also provide constraints on A far more stringent than any current ones, and will 
provide a unique probe of the inflationary observables. Information on the baryon density and Hubble constant will 
complement and perhaps even improve upon current observations. Furthermore, although we have yet to include 
polarization maps in our error estimates, it is likely that they will provide additional information, at least regarding 
ionization history. 

We have attempted to display our results in a way that will be useful for future CMB experimental design. Although 
a satellite mission offers the most promising prospect for making a high-resolution CMB map, our error estimates 
should also be applicable to complementary balloon-borne or ground-based experiments which map a limited region 
of the sky. The estimates presented here can also be used for a combination of complementary experiments. 

Although we have been able to estimate the precision with which CMB temperature maps will be able to determine 
cosmological parameters, there is still much theoretical work that needs to be done before such an analysis can 
realistically be carried out. To maximize the likelihood in a multidimensional parameter space, repeated evaluation 
of the CMB spectrum for a broad range of model parameters is needed. Therefore, quick and accurate calculations 
of the CMB anisotropy spectrum will be crucial for the data analysis. Several independent numerical calculations of 
the CMB anisotropy spectrum now agree to roughly 1% j5^. However, these calculations typically require hours of 
workstation time per spectrum and are therefore unsuitable for fitting data. We have begun to extend recent analytic 
approximations to the CMB spectra ^0 52 1 with the aim of creating a highly accurate and efficient power spectrum 
code. Our current code evaluates spectra in a matter of seconds on a workstation, though our calculations are not yet 
as accurate as the full numerical computations, except in a limited region of parameter space. It is likely, however, 
that the analytic results can be generalized with sufficient accuracy. 

The other necessary ingredient will be an efficient and reliable likelihood-maximization routine. Preliminary fits to 
simulated data with a fairly simple likelihood-maximization algorithm suggest that the cosmological parameters can 
indeed be reproduced with the precision estimated here 

In summary, our calculations indicate that CMB temperature maps with good angular resolution can provide an 
unprecedented amount of quantitative information on cosmological parameters. These maps will also inform us about 
the origin of structure in the Universe and test ideas about the earliest Universe. We hope that these results provide 
additional impetus for experimental efforts in this direction. 
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FIG. 1. Predicted multipole moments for standard CDM and variants. The heavy curves in each graph are for a model with 
primordial adiabatic perturbations with Q, = 1, A = 0, ns = 1, Qth^ = 0.01, h = 0.5, a = 0, and no tensor modes. The graphs 
show the effects of varying Q, A, h, Treion = 0, and Qbh^ while holding all other parameters fixed. In the Q panel, from left 
to right, the solid curves are for Q, — 1, Q = 0.5, and Q — 0.3. The curves in the flth^ panel are (from lower to upper) for 
flbh^ — 0.01, Qbh^ = 0.03, and 0,bh^ = 0.05. In the h panel, the heavy curves is for h = 0.5, while the other two curves are for 
h — 0.3 (the upper light curve) and h = 0.7 (the lower light curve). The curves in the A panel are for (from lower to upper) 
A = 0, A = 0.3, and A = 0.7. 
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FIG. 2. Simulated data that might be obtained with a CMB 
noise level of w'^ = 2 x 10~^®. 
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FIG. 3. The standard errors for Q, A, Qbh^, and h that can be obtained with a full-sky mapping experiment as a function of 

the beam width ^fwhm for noise levels — 2 x 10"^^, 9 x 10^^'', and 4 x 10^^'' (from lower to upper curves). The underlying 
model is "standard CDM." The solid curves are the sensitivities attainable with no prior assumptions about the values of any 
of the other cosmological parameters. The dotted curves are the sensitivities that would be attainable assuming that all other 
cosmological parameters, except the normalization (Q), were fixed. The results for a mapping experiment which covers only a 
fraction /sky of the sky can be obtained by scaling by f^y^'^ [cf., Eq. (20)]. 
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FIG. 5. The standard errors on the inflationary observablos, ns, ut, r = Qt/Q%i and Q, that can be obtained with a full-sky 
mapping experiment as a function of the beam width ^twhm for noise levels = 2 X 10"^^, 9 x IQ-^^, and 4 x lO"^* (from 
lower to upper curves). The parameters of the underlying model are the "standard=CDM" values, except we have set r = 0.28, 
ns = 0.94, and nr = —0.04. The solid curves are the sensitivities attainable with no prior eissumptions about the values of any 
of the other cosmological parameters. The dotted curves are the standard errors that would be attainable by fitting to only 
these four inflationary observables and assuming all other cosmological parameters are known. (Note that this differs from the 
dotted curves in Fig. 4.) The results for a mapping experiment which covers only a fraction /sky of the sky can be obtained by 
scaling by f;^J^ [c.f., Eq. (20)]. 
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FIG. 6. Plots of the log-likelihood £is a function of f2, A, ilbh^, and ns for the "standard-CDM" model with tensor modes 
(so ns = 0.94 and r = 0.28). Note that the log-likelihood looks parabolic. 
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